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TO ALL WHOM IT MAY CONCERN: 

Be it known that, Seemeen S. Karimi, citizen of India and residing in Brookline, MA, and 
Carl R. Crawford, citizen of the United States of America and residing in Brookline, MA, 
have invented certain improvements in a HELICAL INTERPOLATION FOR AN 
ASYMMETRIC MULTI-SLICE SCANNER of which the following description in 
connection with the accompanying drawings is a specification, like reference characters 
on the drawings indicating like parts in the several figures. 
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HELICAL INTERPOLATION FOR AN ASYMMETRIC MULTI-SLICE 

SCANNER 

CROSS REFERENCES TO RELATED APPLICATIONS 

[0001] This application claims the benefit of priority under 35 U.S.C. §1 19(e) from co- 
pending, commonly owned U.S. provisional patent application serial number 60/410,244 
(entitled "Helical Interpolation For An Asymmetric Multi-Slice Scanner" and having an 
attorney docket number ANAK-240PR), which was filed on September 12, 2002. 

FIELD OF THE INVENTION 

[0002] The present invention relates to CT image reconstruction. In particular, the 
present invention relates to helical interpolation for CT image reconstruction from 
asymmetric, multi-slice CT scans. 

BACKGROUND 

[0003] CT scans allow an image of the internal structure of a target object to be 
generated, one cross-sectional slice at a time. Typically, the target object is an 
anatomical region of a patient, although CT systems can also be used in non-medical 
applications, for example explosive detection. In a CT system, x-rays emitted from an x- 
ray source are passed through a region of the object, then are detected by a detector 
assembly. The detector assembly, consisting of one or more rows of detector elements, 
generates detection signals indicative of the attenuated intensities of the x-rays that have 
traversed the object. The detection signals are sent to a computer, which implements 
signal and image processing techniques to reconstruct a tomographic image of the object. 

[0004] In a helical CT scanner, the patient is translated (typically at a constant speed), 
while the x-ray source and the detector assembly rotate around the patient. As the patient 
is moving, the data for the prescribed number of axial slices of the target region (within 
the patient) is acquired. Because the patient table is translated at a constant speed along 
the axis of rotation of the gantry during helical scanning, the image location, also called 
the slice plane, constantly moves in the axial direction. The trajectory of the source 
relative to the slice plane maps out a helix, generating projection data from which axial 
image slices may be reconstructed. Helical scanners offer a number of advantages, 
including reduced scanning time, improved image quality, and better control of contrast. 
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[0005] In a multi-slice helical scanner, more than one row of detector elements are 
arranged side-by-side along the rotation axis of the scanner. For a given view angle, it is 
thus possible to obtain projections that measure attenuation in multiple sections of the 
scanned object in a single sampling period. Multi-slice scanners are becoming the norm 
for medical CT applications. Multi-slice scanners make rapid acquisition of volumetric 
data possible, because of the larger coverage that they provide (compared to the coverage 
provided by single slice scanners), coupled with helical scanning. However, since in 
helical multi-slice scanners the patient is translated for each sample, the projections 
measure attenuation at different positions within the patient. This necessitates extracting 
consistent data sets for each cross-sectional slice position from the helical multi-slice 
data, in order to reconstruct each cross-sectional slice, thereby adding further 
complications to the image reconstruction process. 

[0006] Currently, volumetric reconstruction for multi-slice scanners is typically 
performed by helical interpolation, followed by 2D filtered back-projection. Generally, 
the number of rows in multi-slice scanners is limited to about four, in order to prevent 
image artifacts. When the number of rows is limited to four, it has been found that 
helical interpolation ignoring the divergence of the x-ray beam, coupled with 2-D 
backprojection, is sufficient to provide clinically acceptable image quality. In this case, 
an interpolation filter is used to estimate virtual fan beam data for an image at a desired 
plane, given the positions of the plane, the source and the detector. For a larger number 
of rows, it has been observed that 3-D backprojection is necessary to provide acceptable 
image quality. Asymmetric fan beams have been used to increase a scanner's field of 
view (FOV) in a cost effective manner, because the FOV can be increased by increasing 
the number of detectors on only one side of the fan beam. When the FOV is increased in 
this manner, resolution is lowered in the extended part of the FOV. 

[0007] All of the features described above, namely helical scanning, multiple detector 
rows and asymmetric fan beam, are desirable features for a CT scanner (or PET scanner), 
because of the advantages described above. In order to implement such an asymmetric, 
multi-slice, helical scanner, an appropriate image reconstruction method and system is 
required. 



BST99 1369002-1.056230.0598 



3 



[0008] Accordingly, it is an object of the present invention to provide a method and 
system for CT image reconstruction, that can be used with an asymmetric multi-slice 
helical scanner, and that produces acceptable image quality. 
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SUMMARY 

[0009] The present invention provides a method and system for image reconstruction 
from asymmetric, multi-slice helical scans. In overview, it is proposed that 
complementary interpolation be performed for projection data generated within the 
symmetric region of the x-ray beam in the multi-slice scanner, and that direct 
interpolation be performed for projection data generated within the asymmetric region of 
the x-ray beam. It is also proposed that blending regions be created near the boundaries 
between the symmetric and asymmetric regions, and that a combination of direct and 
complementary interpolation be performed in these blending regions. 

[0010] A CT system in accordance with one embodiment of the present invention 
includes an x-ray source for generating x-rays. The x-ray source is mounted on a gantry 
for rotation about a rotation axis. An x-ray detector system is mounted opposite the x-ray 
source for providing a set of projection data with respect to the object. The CT system 
further includes an interpolator for interpolating the proj ection data from the detector 
array onto a slice plane by multiplying the data with helical interpolation weights. The 
interpolator includes an interpolation weight generator for generating the helical 
interpolation weights. The CT system also includes an image reconstructor for 
reconstructing a tomographic image of the object using the helically interpolated 
projection data. 

[001 1] The reconstructed images are perpendicular to the rotation axis, and therefore 
generally referred to as axial images or axial slices. The plane defined by an axial image 
is generally referred to as the slice plane. A plane parallel to the plane of gantry rotation 
is generally referred to as an axial plane. 

[0012] The x-ray detector system includes a plurality of rows of detector elements. For 
each detector element, a corresponding detector ray is defined by x-ray photons traveling 
from the x-ray source to the detector element. The detector rays in a multi-slice scanner 
define an x-ray beam that is generally referred to as a cone beam. The cone beam of an 
asymmetric multi-slice scanner is asymmetric in an axial plane. The x-ray beam includes 
a symmetric region in which one or more complementary rays can be found for each 
detector ray, and an asymmetric region in which no complementary ray can be found for 
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any detector ray. For a given detector ray, a complementary ray is defined in the present 
application as a detector ray whose projection on an axial plane is anti-parallel to the 
projection of the given detector ray on an axial plane, with the source and detector array 
positions reversed. Complementary rays as defined in the present application are not 
necessarily on parallel lines, but their projections on an axial plane substantially coincide. 

[0013] A set of virtual fan beam data per slice are estimated from the helical cone beam 
data. For projection data resulting from detector rays lying within the symmetric region 
of the x-ray beam, the helical interpolation weights are called complementary 
interpolation weights. The complementary interpolation weights weigh complementary 
projection data from the different detector rows in proportion to the distance from the 
rows to the slice plane. 

[0014] For projection data resulting from detector rays that lie within the asymmetric 
region of the x-ray beam, the helical interpolation weights are called direct interpolation 
weights. The direct interpolation weights weigh the projection data from different rows 
in proportion to the distance from each detector row to the slice plane. 

[0015] For projection data resulting from detector rays that lie within the symmetric 
region of the x-ray beam but are adjacent to the asymmetric region of the x-ray beam, the 
helical interpolation weights are a combination of direct and complementary interpolation 
weights. 

[0016] A method of reconstructing at least one image of an object, in accordance with 
one embodiment of the present invention, includes helically scanning the object with x- 
rays to acquire tomographic projection data representative of the object while the object 
is translated along an axis. The x-rays are generated by an x-ray source mounted on a 
gantry for rotation about the axis along which the object is translated. The x-rays are 
incident upon a multi-slice x-ray detector system having a plurality of substantially 
parallel rows of detector elements. For each detector element, a corresponding detector 
ray is defined by x-ray photons traveling from the x-ray source to the detector element. 
The detector rays define an x-ray beam that is asymmetric in the plane perpendicular to 
the axis of gantry rotation. The x-ray beam includes a symmetric region in which one or 
more complementary rays can be found for each ray, and an asymmetric region in which 
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no complementary ray can be found for any ray. 

[0017] The method further includes helically interpolating the projection data by 
multiplying the data with helical interpolation weights. For projection data resulting 
from detector rays lying within the symmetric region of the x-ray beam, the helical 
interpolation weights are complementary interpolation weights. For projection data 
resulting from detector rays that lie within the asymmetric region of the x-ray beam, the 
helical interpolation weights are direct interpolation weights. For projection data 
resulting from detector rays that lie within the symmetric region of the x-ray beam but are 
adjacent to the asymmetric region of the beam, the helical interpolation weights are a 
combination of direct and complementary interpolation weights. 

[0018] The method further includes reconstructing a tomographic image of the object, 
using the helically interpolated projection data. In one embodiment, a 2D filtered 
backprojection of the helically interpolated projection data may be performed. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0019] FIG. 1 illustrates a schematic block diagram of an asymmetric multi-slice CT 
scanner constructed in accordance with one embodiment of the present invention. 

[0020] FIG. 2 illustrates an asymmetric x-ray beam, used in one embodiment of the 
present invention. 

[0021] FIG. 3 illustrates a helical asymmetric interpolation filter, in accordance with one 
embodiment of the present invention. 

[0022] FIG. 4 illustrates a reconstructed CT image of a cylindrical phantom. 

[0023] FIG.s 5 A, 5B, and 5C illustrate reconstructed CT images of a phantom consisting 
of two spheres. 
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DESCRIPTION 

[0024] The present invention is directed to a helical interpolation filter for reconstructing 
a tomographic image of an object, using an asymmetric, multi-slice helical CT scanner. 

[0025] FIG. 1 illustrates a schematic block diagram of an asymmetric, multi-slice CT 
scanning system 10, constructed in accordance with one embodiment of the present 
invention. The CT scanning system 10 includes an x-ray source 12 mounted on a gantry 
13 for rotation about a rotation axis 15, and an x-ray detector system 17 mounted 
opposite the x-ray source 12 for providing projection data with respect to an object 18, as 
the object 18 is translated along the rotation axis 15. The center of the circle formed by 
the rotation of the x-ray source 12 is the isocenter 1 1 of the CT system 10. The x-ray 
source 12 can be considered as substantially a point source. The detector system 17 
includes a plurality of rows 19 of detector elements. In a multi-slice scanner, such as the 
CT system 10, data from multiple detector rows are used to reconstruct volume images. 

[0026] For each detector element, a corresponding detector ray (or x-ray path) 21 can be 
defined by the x-ray photons traveling from the x-ray source to the detector element. 
Each detector element generates detection signals indicative of the intensity of its 
corresponding detector ray 21. The center ray 25 is the line from the x-ray source 12 
through the isocenter 1 1 . 

[0027] Each individual row 19 of detector elements is typically configured in the shape 
of an arc of a circle. Preferably, the plurality of rows 19 are substantially parallel to each 
other, and are disposed side-by-side along the rotation axis. In this way, during a single 
sampling period, projection data can be acquired that are representative of a plurality of 
sections of the object 18. 

[0028] Because each detector ray 21 is at least in part attenuated by the object it 
encounters in its path, the detection signal from each detector element is representative of 
the attenuation of the portion of the object that lies in the path of the detector ray. The 
raw detection signals are processed to generate a set of projection data, representative of 
the logarithmic attenuation effected by the mass lying in the corresponding detector ray 
path. 

[0029] The CT system includes an interpolator 22 for interpolating the projection data 
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from the detector array onto a slice plane, by multiplying the data with helical 
interpolation weights. The interpolation weights are generated by an interpolation weight 
generator 24, described in detail below in connection with FIG. 3. The CT system further 
includes an image reconstructor 26 for reconstructing a tomographic image of the object, 
using the helically interpolated projection data. In a preferred embodiment of the 
invention, the image reconstructor 26 includes means for performing 2D backprojection 
of the helically interpolated projection data. It should be understood, however, that 
image reconstruction techniques other than 2D backprojection are within the scope of the 
present invention, and that in alternative embodiments of the present invention, the image 
reconstructor 26 may include alternate means for reconstructing images from the 
helically interpolated projection data. 

[0030] Since reconstructed images are perpendicular to the axis of rotation , they are 
generally referred to as axial images or axial slices. In the present application, following 
customary convention, the plane defined by an axial image is referred to as the slice 
plane. A plane parallel to the plane of gantry rotation, i.e. to the plane defined by the 
rotation of the gantry, is referred to as an axial plane. 

[003 1 ] Since the detector rays originate from the point x-ray source and extend at 
different angles to the different detector elements, the multiple detector rays resemble a 
cone. The collection of all the detector rays at a given instant of time is often referred to 
as a "cone beam." If there were a single row of detectors, then the collection of detector 
rays would form a "fan beam," i.e. a fan- shaped beam. The x-ray cone beam can be 
considered to be a collection of fan beams. In the present application, for convenience, 
the term "x-ray cone beam" will be used to refer to the x-ray beam defined by the 
collection of all the detector rays at any point in time, or at any time interval. In the 
present invention, the CT system is an asymmetric CT system, in which the x-ray cone 
beam is asymmetric in an axial plane. 

[0032] Asymmetric fan beams can be used to increase a scanner's field-of-view (FOV) 
with lower image quality at the outer part of the FOV, by increasing the number of 
detectors only on one side of the fan beam. One example of a CT scanner using an 
asymmetric fan beam is the "A" scanner, manufactured by Analogic Corporation. In 
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such scanners, one side of the fan beam covers the entire FOV, while the other side fails 
to cover a portion of the FOV. In the part of the FOV that is covered by only one side of 
the fan beam, the sampling frequency is reduced. The image quality is therefore 
compromised in this part of the FOV. 

[0033] A set of virtual fan beam data per slice are estimated from the x-ray cone beam 
data. For each ray in the virtual fan beam, rays at the same azimuthal angle as the given 
fan beam ray, and complementary rays, are identified in the x-ray cone beam. The 
azimuthal angle is defined in this application as the angle made by the projection of a ray 
onto the plane of gantry rotation with a fixed axis in said plane. The interpolator 
interpolates projection data from different rows at these rays, to estimate the virtual fan 
beam data for the slice. If more than one rotation of views is used to generate one slice, 
then multiple projection rays and multiple complementary projection rays for a given 
azimuthal angle may be obtained and used. 

[0034] A number of reconstruction algorithms are known for multi-slice scanners. Some 
algorithms are able to perform exact reconstruction, while others achieve various degrees 
of approximate reconstruction. Exact reconstruction techniques use 3D backprojection, 
while approximate algorithms may use either 3D or 2D backprojection. When the 
number of rows is limited to four, it has been found that helical interpolation ignoring the 
cone beam divergence, coupled with 2D backprojection, is sufficient to provide clinically 
acceptable image quality. 2D reconstruction assumes that the x-rays incident on all rows 
are perpendicular to the axis of rotation, ignoring the fact that the cone beam is in reality 
divergent. The interpolation filter 22 is used to estimate the data at the slice plane, given 
the positions of the plane, the source and detector. 

[0035] FIG. 2 illustrates a 2-dimensional view of an asymmetric x-ray beam 100, in an 
exemplary asymmetric CT scanner described in conjunction with FIG. 1 . The 2- 
dimensional view is parallel to an axial plane, and therefore appears as a fan-shaped 
beam. The x-ray beam 100 includes a symmetric region 102, in which at least one 
complementary ray can be found for each detector ray within the symmetric region, and 
an asymmetric region 104, in which no complementary ray can be found for any detector 
ray within the asymmetric region. 
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[0036] For a given detector ray, a complementary ray is defined in the present application 
as a detector ray whose projection on an axial plane is anti-parallel to the projection of 
the given detector ray on an axial plane, with the source and detector array positions 
reversed. Complementary rays as defined in the present application are not necessarily 
on parallel lines, but their projections on an axial plane substantially coincide. 

[0037] As illustrated in FIG. 2, there are no detector elements in region 106, thereby 
excluding the possibility of complementary detector ray pairs. Typically, for all detector 
rays lying with the symmetric region 102 of the x-ray beam, the projection data used in 
interpolation for complementary ray pairs may come from different ones of the multiple 
rows of detector elements, or from all the rows. 

[0038] The azimuthal angle of the x-ray source 1 12 position is denoted as p (the "fan 
angle"). The fan angle is defined as the angle between a line from the x-ray source 
(indicated in FIG. 2 by the reference numeral 1 12) to the isocenter of the CT system, and 
a fixed coordinate axis in the gantry rotation plane. The angle made by the projection of 
a detector ray onto an axial plane with the center ray (illustrated in FIG. 1) of the x-ray 
beam is denoted y, and defined as the ray angle in this application. The value of y at the 
boundary of the symmetric region of the x-ray beam is indicated as y s . The value of y at 
the end of the asymmetric region of the x-ray beam is indicated as y a . 

[0039] Helical interpolation is implemented by creating a set of weights, and multiplying 
the weights by the projection data. FIG. 3 schematically illustrates a helical asymmetric 
interpolation filter, in accordance with one embodiment of the present invention. In the 
illustrated embodiment, the helical interpolation filter is designed for a four-row scanner 
with an asymmetric detector array; however the method and system of the present 
invention may be used for CT scanners having a number n of rows, where n is different 
from four. 

[0040] FIG. 3 schematically illustrates an symmetric region 200, an asymmetric region 
202, a region 203 in which no data are collected, and a pair of blending regions 204, 
within the asymmetric x-ray cone beam data. The symmetric region of the x-ray cone 
beam comprises a range {-y s , y s ) of ray angles, as seen in FIG. 3. The symbol y b is the 
value of y at the start of the blending region within the symmetric region of the x-ray 
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cone beam. The fan angles fi t indicate the view angle at which the row / is directly under 
the slice plane at isocenter. The symbol p* denotes the ray complementary to the angle 
fii. It depends upon the value of y for each ray. The superscript indicates whether the 
complementary ray comes from the current rotation (+) or the previous rotation (-). The 
interpolation weight for each row x is denoted as w,. A weight generated by 
complementary interpolation is labeled w" and a weight generated by direct interpolation 
is labeled wf . 

[0041] In the helical interpolation filter of the present invention, complementary 
interpolation is performed for data that lie within the symmetric region 200 of the x-ray 
cone beam, because in this region all rays have at least one complementary ray that 
results from reversing the positions of the source and the detector array. In the 

0 

asymmetric region 202, there are no complementary rays. In this region, direct data 
interpolation is used. 

[0042] A blending of interpolation weights across the direct and complementary 
interpolation portions is necessary to ensure continuity of the weights across the 
boundaries between the symmetric and asymmetric regions. Therefore, a blending region 
204 is created near the boundary between the complementary and direct interpolation 
regions. In the blending regions 204, a combination of the direct and complementary 
interpolation is used. Since there are no complementary ray pairs in the asymmetric 
region 202 of the detector, blending is performed within the symmetric region 200. 

[0043] For complementary interpolation, complementary ray pairs are identified for rays 
in the virtual fan beam. Complementary rays are related by the following equation. 

Pc = P + tc + 2y 
7c = -Y 

where p c indicates the fan angle of the complementary ray, and y c indicates its ray angle. 

[0044] In general, each sample of the complementary data pair comes from different 
rows. The complementary interpolation weights weigh complementary projection data 
from different detector rows in proportion to the distance from the rows to the slice plane. 
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Complementary rays are not truly identical, due to cone beam divergence in the axial 
direction. Further, due to quarter detector offset, the projection of complementary rays 
onto the axial plane do not truly coincide. Therefore, in the implementation of helical 
interpolation, the projection data are weighted according to their positions in the x-ray 
cone beam, but not combined. 

[0045] In the illustrated preferred embodiment, the complementary interpolation weights 
are calculated according to the following equations, using the conventions described in 
paragraph 40: 
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[0046] The complementary weights should be independent of the absolute view angles. 
By way of example, in a scanner mode in which p, = , the /?- values for the other rows 

(i = 2, . . . 4) are given as follows: p 2 = -f- , P 3 = In , and P 4 = -y 1 . These values of 



P,. (i= 1, . . . , 4) should be used to evaluate the weight formulas above. In this case, the 
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range of the p variable is from 0 to 



IOjt 



[0047] For the asymmetric region 202 of the x-ray cone beam, direct interpolation is 
performed. The angular range in which direct interpolation is performed is given by 
(Y^'Ya) ' where y b is the ray angle in the symmetric region at which blending is started, 

as shown in Figure 2. The symbol y a is the maximum ray angle of the longer side, as 

shown in Figures 2 and 3. 

[0048] Direct data interpolation means that for a ray in the virtual fan with a given source 
fan angle and a given ray angle, the data from at least two different rows in the x-ray 
cone beam are interpolated. Direct interpolation weights weigh the data from different 
rows at a given source fan angle and ray angle, in proportion to the distance from each 
detector element row to the slice plane. The disadvantage of direct data interpolation is 
that the sampling frequency in the axial direction is reduced and the interpolation 
distances are larger. This means that the slice sensitivity profile (SSP) is degraded. 

[0049] In the illustrated embodiment, in which the number of detector element rows is 
four, the direct interpolation weights are defined by the following relationship, using the 
conventions described in paragraph 40: 
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In equations (6) - (9) above, a, x\ 9 andx 4 are defined as in equation (5). 



[0050] As seen in FIG. 3, the blending regions 204 are defined within the symmetric 
region 200, one blending region being located adjacent the boundary between the 
symmetric region 200 and the asymmetric regions 202, the other blending region being 
located adjacent the boundary between the symmetric region 200 and the region 203 in 
which no data are collected. Blending smoothes the transition between the asymmetric 
and symmetric regions. Blending weights are applied within the symmetric part of the 
boundary, so the complementary data must also be weighted accordingly. The blending 
regions 204 extend between angle ranges (y b , y s ), and (-y b ,-Ys) ? where Yb represents the 
value of y at the start of the blending region within the symmetric region 200 of the x-ray 
beam, and y s represents the value of y at the boundary between the symmetric region 200 
and the asymmetric region 202. 

[0051] For projection data resulting from rays that lie within the blending regions 204, 
the interpolation weights are blending interpolation weights, which are a combination of 
direct and complementary interpolation weights. In the illustrated embodiment, the 
blending interpolation weights are defined by the following relationship: 
within the range (y a> y f ), 

w, (P ,Y ) = (1 - a t WK(P ,y ) + a f (x)wf (P ,y ) ; (10) 
within the range (^y s ,-y b ) , 



H 4 (P > y) = a-a / WKO,y); 



(11) 
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Y s-Jb 

[0052] In helical interpolation, the angular range of the data that is weighted is greater 
than 2n. The redundant data are added so that the output data spans an angular range of 
27i. Rebinning the set of 2n data fills in the missing data on one side of the fan beam. 
The rebinning process performs tangential interpolation, providing hybrid views. Each 
hybrid view contains rays that are parallel in the angular direction, but are not equally 
spaced. Each hybrid view is added to its opposite view, located n away. The opposite 
view is flipped before adding. The asymmetric data from the opposite view thus fill in 
the missing data in each view. 

[0053] FIG. 4 illustrates a reconstructed CT image of a simulated cylinder phantom, used 
to test the helical interpolation filter described in connection with FIG. 3. The phantom is 
a cylinder whose axis is parallel to the axis of rotation. The phantom is placed at the 
boundary between the symmetric and asymmetric regions. For comparison, an image 
using a symmetric fan beam that covers the full field of view is also shown in FIG. 4. 

[0054] The cylinder phantom is substantially uniform in the axial direction, so it allows 
us to verify whether the helical interpolation weights for each ray are correctly 
normalized and are continuous. Figure 4 shows a substantially uniform attenuation 
within the cylinder. There are no artifacts at the boundary of the symmetric and 
asymmetric regions, which means that the transition is continuous, and both regions show 
the same CT value, which means that the normalization of weights is correct. The 
interpolation weights for each ray add to one. 

[0055] FIG.s 5 A, 5B, and 5C illustrate reconstructed CT images of another phantom, 
used to test image quality. The phantom in each of the figures consists of two spheres of 
radius 30 mm. In these figures, the window width is 50 HU (Hounsfield units). The 
sphere at the center of the image is centered on the axis of rotation, and the other sphere 
in each image is centered at the boundary of the asymmetric and symmetric regions. 
FIG. 5 A is the image produced using an asymmetric x-ray fan beam. For comparison, an 
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image using a symmetric fan beam that covers the full field of view is also shown. FIG. 
5B is the image produced using a symmetric x-ray fan beam. FIG. 5C shows the 

difference image. As seen from FIG.s 5A - 5C, some artifacts are observed. 

« 

[0056] The artifacts in the background (as seen in FIG.4 and in FIG.s 5A - 5C) are 
caused by linear helical interpolation and also because the cone beam was ignored. The 
artifacts in the symmetric fan image are slightly smaller because complementary data 
extends to the extreme ends of the beam. The linear interpolation distances are smaller 
for the outer detector elements in the fully symmetric beam, as compared to the 
interpolation distances for the detector elements in the asymmetric beam. However, the 
artifacts from the off-center object are larger than they are from the centered object 
because the interpolation distances are larger than they are at isocenter, as seen in Figure 
3. At the maximum ray angle, the interpolation distance equals one, which is the same 
distance as in direct data interpolation. 

[0057] In sum, the present invention provides a method and system for reconstructing 
helical, multi-slice data from an asymmetric beam. The images created using the 
technique of the present invention have a satisfactory image quality when tested using 
phantoms. 

[0058] While the invention has been particularly shown and described with reference to 
specific preferred embodiments, it should be understood by those skilled in the art that 
various changes in form and detail may be made therein without departing from the spirit 
and scope of the invention as defined by the appended claims. 
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